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SUMMARY 

A segment of the heat exchanger proposed for use the In the NASA Lewis 
Altitude Wind Tunnel (AWT) facility has been tested under dry and Icing condi- 
tions. The heat exchanger has the largest pressure drop of any component In 
the AWT loop. It Is therefore critical that Its performance be known at all 
conditions before the final design of the AWT Is complete. 

The heat exchanger segment was tested In the NASA Lewis Icing Research 
Tunnel (IRT) In order to provide an Icing cloud environment similar to what 
will be encountered In the AWT. 

Dry heat transfer and pressure drop data were obtained and compared to 
correlations available In the literature. The effects of Icing sprays on heat 
transfer and pressure drop were also Investigated. 


INTRODUCTION 

A schematic of the NASA Lewis Altitude Wind Tunnel (AWT) and proposed 
modifications Is shown In figure 1. The heat exchanger must remove the heat 
added to the air stream by the fan and other systems necessary to conduct tests 
of engines. Inlets, props, and airframes at altitude and In Icing conditions. 

A preliminary engineering report (PER) has been written by the Sverdrup Corpo- 
ration under contract NAS3-24024-AE. A heat exchanger configuration has been 
recommended by the contractor. 

The heat exchanger has the largest pressure drop of any component In the 
AWT loop. It Is therefore critical that the performance of the heat exchanger 
be known at all operating conditions before the final design of AWT Is 
complete. 

The purpose of this report Is to present results of heat transfer and 
pressure drop tests conducted In the NASA Lewis Icing Research Tunnel on a 
segment of the recommended heat exchanger. The effects of Icing sprays on 
heat exchanger performance were also Investigated. 


SYMBOLS 

Af heat exchanger frontal area (cross-sectional area of duct upstream of 

heat exchanger) 




Af in 
Ahx/Ats 

Am 

A 0 

Ar 

d 

F 

ha 

h R 

^alr 

LMTD 

LWC 

m 

Pr 

Ps 

Ptl 

Q 

q 

R 

Re 

r h 

fblk 

Ti 

Tr 

u 

X 


total fin heat transfer area 

ratio of duct cross-sectional area upstream of heat 
exchanger to tunnel test section cross-sectional area 

minimum free flow area in heat exchanger 

total air side heat transfer surface area 

refrigerant side heat transfer area 

diameter of heat exchanger tube 

least squares function defined by equation (12) 

air side heat transfer coefficient 

refrigerant side heat transfer coefficient 

thermal conductivity of air 

log mean temperature difference defined by equation (5) 

liquid water content 

mass flow rate In heat exchanger duct 

Prandtl number of air 

average static pressure 

total pressure at rake tube 1 

heat exchanger total heat load defined by equation (4) 

dynamic pressure upstream of heat exchanger defined by equation (9) 

gas constant for air 

Reynolds number based on heat exchanger tube diameter defined by 
equation (6) 

hydraulic radius of heat exchanger air side 

bulk temperature of air upstream or downstream of heat exchanger 
temperature nearest rake total pressure tube 1 
refrigerant temperature 

overall heat transfer coefficient defined by equation (7) 

fraction of LWC removed upstream of the heat exchanger by turning 
vanes, fan, etc. 
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velocity calculated for rake total pressure tube 1 
Ap pressure drop across heat exchanger 

n air side surface efficiency defined by equation (14) 

nfin fin efficiency 

v viscosity 

p \ density at rake total pressure tube 1 

Subscripts 


AWT 

IR1 

SEG 

1 

2 


refers to the Altitude Wind Tunnel 
refers to Icing Research Tunnel 

refers to heat exchanger segment tested In Icing Research Tunnel 
refers to conditions upstream of the heat exchanger 
refers to conditions downstream of heat exchanger 


AWT REQUIREMENTS 
Heat Transfer 

AWT heat exchanger requirements were taken from the 50 percent PER- 
sectlon 9 process systems (ref. 1). Heat loads Include the fan, steam heated 
turning vanes, and the heat load due to Ice formation. Table I shows the heat 
loads, flow rates, air temperatures, Reynolds numbers, and the required heat 
exchanger area-overall heat transfer coefficient product. The area-overall 
heat transfer coefficient product was calculated assuming a 225 K (-55 °F) 
constant temperature flooded coll heat exchanger. Reynolds number was calcu- 
lated based on free stream conditions and a 2.54 cm (1.0 In) diameter tube. 

The maximum required overall heat transfer coefficient occurs at the 
0.7 Mach number, 10 058 m (33 000 ft) altitude, 244 K (-20 °F) condition. The 
greatest heat load does not occur at this condition, but the combination of 
heat load and temperatures required make this the most demanding point for 
heat exchanger heat transfer performance. 


Pressure Orop 

The maximum heat exchanger pressure drop occurs at the point of maximum 
flow rate which Is the 0.52 Mach number, 1524 m (5000 ft), 293 K (68 °F) 
condition. 
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Icing 


The test section Icing environment goals for the AWT are: 

Static temperature, K (°F) 233 to 273 (-40 to +32) 

3 

Liquid water content, gm/m . .. . 0.1 to 3.0 

Volume mean droplet size, ym 5 to 50 

Figure 2 shows the present FAA FAR 25 Icing envelope along with the AWT 

goal . 


APPARATUS 

Icing Research Tunnel 

The heat exchanger segment was tested In the NASA Lewis Icing Research 
Tunnel (IRT). This facility was chosen because It could provide the necessary 
cold temperatures and Icing sprays. A circuit diagram of the IRT with the 
location of the heat exchanger segment tested Is shown In figure 3. This 
location was chosen because water sprayed from the nozzles upstream of the 
contraction, must travel through the test section, two sets of turning vanes, 
and the main drive fan before entering the heat exchanger just as It will do 
In the proposed AWT configuration. 


Refrigeration System 

The heat exchanger segment was connected In parallel with the IRT refrig- 
eration system as shown In figure 4. Type E thermocouples and strain gauge 
pressure transducers were Installed In the Inlet and outlet refrigerant lines. 
Refrigerant 12 Is used In the IRT system. The Inlet line was a 2 In copper 
tube and the control valve had a C v of 12. The proposed AWT refrigeration 
system will have a 4:1 overfeed; refrigerant will exit the heat exchanger as a 
mixture of three parts saturated liquid and one part vapor. The refrigerant 
exited the heat exchanger segment tested In IRT with several degrees of super- 
heat Indicating that this valve trim was too small. It was thought that this 
did not seriously affect the overall heat transfer and pressure drop data. 


Heat Exchanger 

The heat exchanger segment tested was a refrigerant 12 - air heat 
exchanger constructed of 2.601 cm (1.024 In) outside diameter copper tubes 
(0.124 cm (0.049 In) wall thickness) with seven helically wound fins per Inch. 
The fins were 0.041 cm (0.016 In) thick and 1.52 cm (0.60 In) high. The heat 
exchanger segment was 91.44 cm (36 In) wide by 89.92 cm (35.4 In) high. The 
tubes were arranged In a triangular array 6 rows deep with 15 tubes per row. 
Spacing between the tubes In each row was 5.748 cm (2.263 In). The heat 
exchanger segment parameters are given In table I. 

Three type E thermocouples were attached to the Inlet side of the 2nd, 
7th, and 13th "U"-tubes (up from the bottom) that connect the 1st and 2nd rows 
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of heat exchanger tubes. Three similar thermocouples were attached to the 
outlet side of the 3rd, 8th, and 14th tubes of the last tube row. The thermo- 
couples were held to the outside of the tubes with hose clamps. 


Heat Exchanger Duct 

The heat exchanger was Installed In the 91.44 x 89.92 cm (36 x 35.4 In) 
plywood duct shown In figure 5. The duct extended 105.4 cm (41.5 In) upstream 
and downstream of the heat exchanger segment. The Inlet of the duct was 
rounded with a circular contour of 15.88 cm (6.25 In) radius. Air was pulled 
past the heat exchanger segment by a centrifugal fan which had axial flow Inlet 
and outlet. The fan had a set of damper vanes at the Inlet so the flow could 
be adjusted from zero to the maximum of about 15.4 kg/sec (34 Ibm/sec). The 
test duct Is shown schematically In figure 6. 

When the system was first checked out an acoustic resonance was found to 
exist with the fan damper set at about the 50 percent position. This was par- 
tially eliminated by Installing a plywood splitter In the duct downstream of 
the heat exchanger. The splitter was horizontal and divided the duct In half. 

The mass flow rate of air In the duct was measured with total pressure 
rakes constructed as shown In figure 7. Each rake contained 16 total pressure 
tubes and 8 type E thermocouples. The position of the total pressure and 
thermocouples tubes corresponded to the centers of equal areas for the rectan- 
gular duct. There were two Identical rakes upstream and two rakes downstream 
of the heat exchanger. On the upstream side of the heat exchanger the rakes 
were Installed so they bisected the duct In the horizontal and vertical direc- 
tions. Downstream of the heat exchanger the vertical rake bisected the duct 
but the horizontal rake was mounted at a point 3/4 of the duct height up from 
the bottom. This was necessary due to the splitter downstream of the heat 
exchanger. When not taking data, the rake total pressure tubes were discon- 
nected from the pressure measuring system and purged with dry air to keep out 
water and Ice. This was accomplished with a shuttle valve; there were not 
enough of these valves available to purge the static taps In a similar manner. 

Thermocouples were all type E and were referenced to a 339 K (150 °F) 

oven. 


Static pressures were measured at both the upstream and downstream rake 
stations. Eight static taps were located at each station. The static taps 
were located ±15.2 cm (6 In) on either side of where the rake penetrated the 
duct wall . 

Both static and total pressures were measured using an electronically 
scanned pressure (ESP) system. Each pressure was measured by a separate 
±0.007 MPa (1.0 psld) transducer which was read and digitized once every sev- 
eral seconds by a minicomputer. Every 350 scans (about every 10 min) the 
transducers were automatically calibrated by mechanically connecting them to 
three separate reference pressures. This was accomplished automatically with 
a shuttle valve. The three reference pressures were atmospheric pressure and 
about 0.003 MPa (0.5 psl) above and below atmospheric pressure. The reference 
pressures were measured with precision digital quartz transducers. 
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All data were recorded on the laboratory data acquisition system which Is 
described In reference 2. For each data point, 20 "scans" or Individual read- 
ings of each data channel were recorded. These 20 scans were then averaged to 
give a single reading for each channel. This helped to eliminate some of the 
noise present In the system. Because of the electronic multiplexing of the 
Individual data channels to the analog to digital converter, there Is a certain 
offset between data channels. In an extreme case this caused a 0.6 °C Indi- 
cated difference In thermocouples that were at the same temperature. This was 
corrected by shorting all the Inputs to the data acquisition system after a 
run, recording this value, then subtracting this "zero" from each reading. 


Data Reduction 

Ouct mass flow rate . - The velocity for each total pressure tube was 
calculated from the pitot tube equation: 



where was the total pressure from rake tube 1, P s was the average 

of the eight static pressure taps, and p\ the density calculated from the 
perfect gas law using the static pressure and the rake temperature nearest 
total pressure tube 1. For example, referlng to figure 6, the temperature 
from thermocouple 5 would be assigned to total pressure tubes 9 and 10. Mois- 
ture In the air due to water sprays was neglected In the calculations. The 
total mass flow rate at each rake was then calculated from 

A 32 

fc.jf E ,,V, (2) 

where Af Is the cross-sectional area of the duct (0.822 m 2 (8.85 ft 2 )). 

Bulk temperature . - The upstream and downstream bulk (or mass averaged) 
temperatures were computed from 


A 32 

Tb,k ' 3iV 1?! PlV ' 

where T^ Is the temperature nearest total pressure tube 1. 

Heat load . - The heat exchanger heat load was calculated from 


(3) 


0 = m i c p ( T blkl “ T 


blk2 J 


( 4 ) 


where irq Is the mass flow rate from the upstream rakes, Cn Is the specific 
heat of air, Tbiki Is the upstream bulk temperature and Tbik2 Is the 
downstream bulk temperature. 


6 



Log mean temperature difference (LMTO) . - The LMTD was calculated from 



where Tr-j and Tr 2 are the Inlet and outlet refrigerant temperatures. 

Reynolds number . - Reynolds number for the heat exchanger was calculated 

from 


Re 


<W d 

y 


( 6 ) 


where m Is the mass flow rate from the upstream rake, d is the diameter of 
the heat exchanger tube (2.54 cm (1 In)), A m Is the minimum free flow area 
In the heat exchanger, and y Is the viscosity evaluated at the upstream bulk 
temperature, Tbiki- 


Overall heat transfer coefficient . - The product of overall heat transfer 
coefficient and total heat transfer area was calculated from 


UA 0 = Q/LMTD (7) 

Pressure drop . - The dimensionless pressure drop was calculated from 

P - P 

AP = ^1 S2 

q q 

P S 1 and P S 2 were the average upstream and downstream static pressures 
respectively, and q was the dynamic pressure upstream of the heat exchanger 
calculated from 



(9) 


Af Is the duct cross-sectional area Immediately upstream of the heat 
exchanger and R Is the Ideal gas constant for air. 


ERROR ANALYSIS 

An error analysis was performed using the method of Kline and McCUntock 
(ref. 3). This method requires the estimation of errors In each measured 
quantity then combines these estimates Into an error estimate for the final 
parameters of Interest, i.e., Reynolds number, Ap/q. etc. Pressure measure- 
ment error for the ESP system was estimated as ±20.7 Pa (±0.003 ps 1 ) based on 
experience with the system. The type E thermocouples were estimates to have a 
calibration error of ±0.25 percent of the difference between the measured 
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temperature and the reference temperature. Errors In the reference tempera- 
ture were not considered Important. The error In calibration of the strain 
gauge transducers used for refrigerant pressure measurements was estimated to 
be ±345 Pa (±0.05 psla). 

In addition to these Instrument errors, there were errors due to the 
digitization and electronic noise In the ESCORT system (ref. 2). 

These errors were estimated for each data channel by calculating the 
standard deviation of the 20 scans for each channel. The total error estimate 
for each data channel was then computed as: 

Total error = V( Instrument error) 2 + (std. dev) 2 (10) 

There Is no error analysis that can account for systematic errors In the 
measurements. In these tests water sprays, frost, and Ice continually contam- 
inated the static pressure taps. Engineering judgement was used to eliminate 
"bad" static tap readings from the data. The total pressure rakes were purged 
with air when not In use and were not subjected to this problem. 

The average errors for all the data points were found to be: 


Reynolds number, percent 7.4 

Overall heat transfer coefficient, percent 11.8 

Ap/q, percent ..... 4.5 


Errors at the low flow rates were much larger (see table III). This was 
because the total pressures were too small to be measured accurately by the 
pressure measuring system. 


RESULTS AND DISCUSSION 

Overall heat transfer and air side pressure drop characteristics of a heat 
exchanger proposed for use In the NASA Lewis Altitude Wind Tunnel (AWT) have 
been tested under dry and Icing conditions In the Icing Research Tunnel. 


Dry Pressure Drop 

The ratios of heat exchanger static and total pressure drops to dynamic 
pressure upstream of the heat exchanger as a function of Reynolds number are 
shown on figure 8. The asterisk and octagon symbols represent data that were 
taken before any Icing sprays and for less than 6 min accumulation of frost due 
to condensation of humidity from the air. The square and diamond symbols rep- 
resent data taken with dry air which allowed longer periods of running with 
virtually no frost buildup on the colls. This was accomplished by cooling the 
IRT to 244 K (-20 °F) and running It at a moderate speed for about 1/2 hr. The 
IRT heat exchanger removed moisture from the air and lowered the dew point to 
around 242 K (-25 °F). Static pressures, represented by the asterisk and 
square symbols, were the most difficult quantities to measure accurately due 
to moisture. Ice, and frost buildup on and In the taps. For this reason the 
total pressure drop was computed by averaging the pressures from the upstream 
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and downstream rakes and subtracting these two averages. These are the octagon 
and diamond symbols. Good agreement Is obtained between total and static 
pressure drops. It Is not known why there Is such a large scatter between the 
data taken with humid and dry air at the higher Reynolds numbers. The data of 
figure 98(a), reference 4 most closely matched the AWT heat exchanger param- 
eters. This data was curve fit and converted to the Reynolds number and pres- 
sure drop-dynamic head ratio used In this report. The resulting relation, 


f = 347.8 Re' 0 ' 265 


( 11 ) 


Is also shown In figure 8. In the Reynolds number range of Interest, the data 
of figure 98(a), reference 4 predicts a much higher pressure drop than was 
measured. Also shown on figure 8 Is a curve obtained from a proprietary cor- 
relation for the pressure drop across a bank of finned tubes. The two corre- 
lations predict nearly the same pressure drop. 


Dry Heat Transfer 

The AWT heat exchanger was designed to have a 4:1 refrigerant overfeed, 
l.e., refrigerant exits the heat exchanger as three parts liquid and one part 
vapor. It was not possible to supply enough refrigerant to the test segment 
to obtain this condition. Refrigerant entered the heat exchanger as a sub- 
cooled liquid, boiled to a vapor, and exited the heat exchanger with several 
degrees of superheat. This was thought to be acceptable for the purposes of 
measuring the dry heat transfer performance and Investigating Icing behavior. 

The multiple refrigerant phases present In the heat exchanger plus the 
lack of a refrigerant flow rate measurement made It Impossible to measure 
directly a refrigerant side heat transfer coefficient. The only heat transfer 
measurements possible were the overall heat transfer coefficient. The overall 
heat transfer coefficient data based on air side area are shown In figure 9. 
The asterisk symbols are the humid air data and the squares represent the dry 
air data. 

To compare the measured data with existing correlations for air side heat 
transfer, the refrigerant side heat transfer coeffflclent must be known. As 
explained previously. It was not possible to measure this directly. A refrig- 
erant side heat transfer coefficient was estimated from the data by minimizing 
the least squares function 


n 

F - H < U meas ' "ca/ (12 > 

1=1 

where U meas are the n measured overall heat transfer coefficients and 
U C al the corresponding calculated overall heat transfer coefficients 
(neglecting tube wall resistance). 


U 


cal 


h a n 


h R A R 


(13) 
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h a Is the air side heat transfer coefficient calculated using a proprietary 
correlation for a bank of finned tubes, hR Is the refrigerant side heat 
transfer coefficient to be determined. A 0 and Ar are the air and 
refrigerant side areas respectively, n is the air side efficiency 

n . [A 0 - A f1n <l-„ f1n )]/A 0 04) 


where Af^ n Is the total fin area and nfi n the fin efficiency deter- 
mined by the method In reference 5. The function, F, Is minimized by varying 
hR. The best estimate for hR was found to be 1883 W/m 2 °C (332 Btu/hr- 
ft 2 -°F). This seems to be a reasonable value considering that It must be an 
average over a tube which contains refrigerant In the liquid, two phase, and 
vapor states. Also shown on figure 9 Is the overall heat transfer coefficient 
calculated using the proprietary correlation and the data In figure 98(a) of 
reference 4 to calculate h a . 


A curve fit of the data In figure 98(a) reference 4 gave 


h 


a 


0.032 


^ Re°' 67 Pr 1/3 
4r h 


(15) 


where k a ^ r Is the thermal conductivity of air evaluated at the mean tem- 
perature, rj, Is the hydraulic radius, and Pr Is the Prandtl number for 
air. hR was fixed at 1883 W/m 2 °C (332 Btu/hr-ft 2 -°F) . There Is good 
agreement between the calculated and measured overall heat transfer 
coefficients. 


Comparison of Dry Measurements with AWT Requirements 

Pressure Drop . - The maximum pressure drop for the AWT heat exchanger will 
occur at the highest mass flow rate. This will be at the 0.52 Mach number at 
1528 m (5014 ft) altitude (table I). At this condition, the Reynolds number 
for the heat exchanger, based on tube diameter and upstream conditions. Is 
37 740. From figure 8 the maximum Ap/q at this Reynolds number would be about 
16. The dynamic pressure, q, would be 


then 


1 / l!806\ 2 53.35(528) 27.68 

q = 2 \3930 / 2116.3 32.2(144) 


0.36 In H 2 0 


(16) 


Ap = 16(0.36) = 5.7 In ^0 (17) 

This pressure drop Is for the heat exchanger only and does not Include the 
additional "blister" losses (l.e., losses due to the extra wide angle diffuser 
upstream of the heat exchanger, see figure 1.). 


Heat transfer . - Maximum required heat exchanger overall heat transfer 
coefficient - area product, UA Is 1.445xl0 6 W/°C (2.7406x106 Btu/hr-°F) 
and occurs at the 0.7 Mach number 33 215 ft altitude condition (table I). The 
Reynolds number at this condition based on tube diameter and conditions 
upstream of the heat exchanger Is 19 615. From the lower curve on figure 9, 
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the overall heat transfer coefficient at this Reynolds number Is 1.308 W/m 2 °C 
(7.43 Btu/hr-ft-°F) . To scale to AWT heat exchanger area the air side area of 
the segment tested In IRT will be multiplied by the ratio of face areas l.e., 

(UA) amt . <UA) SE8 - (7.43)0046.15) 3||2aiil (18 

This relation gave a ( UA ) AWT °f l-82xl0 6 W/°C (3.45x10^ Btu/hr-°F) which 
Is 26 percent larger than the required UA. 


Frost Buildup 

The effect of running humid air over the cold heat exchanger colls was 
Investigated. The heat exchanger colls were kept at 241 K (-25 °F) and air 
flow was maintained at the maximum speed. The dew point of air was measured 
to be 273 K (32 °F). After approximately 30 min of running the colls were 
nearly blocked with frost buildup. This buildup was virtually uniform from 
front to back of the heat exchanger. The cold colls of the heat exchanger 
make a very efficient dehumidifier; It Is therefore recommended that AWT be 
started with dry air, to minimize loss of performance due to frost buildup. 


Icing 

Heat transfer and pressure drop data for the AWT heat exchanger segment 
were taken for several Icing conditions. Three test section total tempera- 
tures, 244, 255, and 264 K (-20, 0, and +15 °F) (nominal temperatures measured 
upstream of the IRT test section) were run as well as two drop sizes, 15 and 
30 ym. All Icing tests were conducted at a 1.36 gm/m 3 LWC In the IRT test 
section. To achieve the same LWC at the actual AWT heat exchanger a higher LWC 
In the AWT test section would be required due to the difference In the heat 
exchanger-test section area ratios between the IRT and AWT as well as the 
anticipated difference In the fraction of cloud removed upstream of the heat 
exchanger between the IRT and AWT. The values of AWT test section LWC that 
would produce the same LWC at the heat exchnager as those tested would be 
1.49 gm/m 3 for the 15 ym drops and 1.71 gm/m 3 for the 30 ym drops. These 
values were estimated using the relationship 


(LWC) awt 


(LWC) 


^ A hx /A TS^AWT 
IRT < A hx /A TS>IRT 


(1 ~ X IRT ) 

^- X AWT> 


(19) 


where (LWC)jrj Is the liquid water content In the IRT test section, (Ah x / 

Ats) AWT and (Ahx/Ais) IRT are the area ratios between the heat exchanger and 
test section cross sections for the AWT and IRT respectively (12.51 and 13.96 
respectively). Xjrj and are the fraction of the Icing cloud 

removed by the turning vanes, screens and fan for the IRT and AWT respectively. 
These were estimated from figure 10(a), (b), which were taken from reference 6. 
These two test points are shown by the symbols on figure 2. Tests were also 
conducted with and without refrigerant flowing through the heat exchanger to 
determine the effect of surface temperature on Ice growth rate. 
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Ice buildup . - For all conditions tested. Ice buildup on the heat 
exchanger colls had the same characteristics. Ice formed only on the leading 
edges of the fins and tubes that were In direct llne-of-slght from upstream. 

Ice formed on the second row of tubes only on the fin leading edges directly 
downstream of the gap between tubes In the first row. The Ice accretion can 
be seen In figure 11. Some of the Ice was knocked off before the photograph 
was taken to give a better Idea of the quantity of Ice accreted. This photo 
was taken after 50 min spray at 264 K (+15 °F) with 15 ym drops at an AWT test 
section equivalent 1.49 gm/m 3 LWC. 

Effect of air temperature . - Total and static pressure drop - dynamic 
pressure ratios are plotted versus Ice accretion time for 244, 255, 264 K 
(-20, 0, and +15 °F) test conditions In figure 12. For all these data the IRT 
was run at a constant test section speed of 67 m/s (150 mph), the drop size was 
15 ym and the liquid water content was 1.36 gm/m 3 In the IRT test section 
(AWT equivalent LWC Is 1.49 gm/m 3 ). Although there Is some scatter In the 0 
and 15 °F data the trend Is clearly toward faster Ice buildup at colder tem- 
peratures. At the 244 K (-20 °F) temperature the pressure drop-dynamic head 
ratio Increases by a factor of 2.3 In 15 min of Icing. During Icing tests the 
fan In the AWT heat exchanger duct was run full speed with the fan damper wide 
open. As Ice built up the flow rate and thus Reynolds number decreased as 
shown on figure 13. 

Effect of drop size . - Figure 14 shows the pressure drop-dynamic head 
ratio for drop sizes of 15 and 30 ym. It Is clear that smaller drop sizes 
cause the heat exchanger to plug with Ice In a shorter time. As seen from 
figure 10(a), the larger the drop size the more water Is removed before reach- 
ing the heat exchanger; thus, the LWC at the heat exchanger will be lower for 
larger size drops. This Is opposite to the trend shown In figure 14. During 
the heat exchanger tests, LWC was measured at the heat exchanger location. 

This Involved measuring the drop size using the soot slide technique. With 
this technique a slide that was coated with lamp black (soot) was exposed for 
a brief time to the air stream containing the water drops. The drops Impinge 
on the slide and cause a "crater" In the soot. The diameter of these craters 
can be correlated to drop size. With the 30 ym drops the craters were round 
Indicating spherical liquid water drops were still present. With the 15 ym 
drops many of the craters were jagged as If Ice crystals had struck the soot 
Instead of liquid water. The presence of all ready frozen Ice crystals could 
be the cause of the formation of less dense Ice and thus faster buildup with 
smaller drop sizes. 

Effect of surface temperature . - The effect of heat exchanger surface 
temperature Is shown on figure 15. Two temperatures were obtained by flowing 
or not flowing refrigerant through the heat exchanger. For the no refrigerant 
case the surface temperature was determined by the air temperature and was 
measured to be about 0.6 °C (1 °F) higher than the air temperature. With 
refrigerant flowing the difference between air temperature at the Inlet to the 
heat exchanger and the first row coll surface temperature was about 6.7 °C 
(12 °F) . A discontinuity Is apparent In the data with refrigerant off. It Is 
believed that the data after 10 min of Icing spray with the refrigerant off Is 
In error and that there Is no effect of heat exchanger surface temperature 
(unless of course the air Is above freezing). This Is substantiated by visual 
observation of the Ice buildup; no difference could be seen between the two 
cases for the same Ice accretion times. 
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Heat transfer . - The effect of Ice accretion on overall heat transfer 
coefficient Is shown on figure 16. This test was conducted at 244 K (-20 °F), 
15 vm drops and an equivalent AWT test' section liquid water content of 
1.49 gm/nw. The Initial Reynolds number for this test was 57 800. This 
corresponds to full flow for the heat exchanger duct. As Ice built up on the 
colls the flow rate and thus Reynolds number decreased (fig. 13). The heat 
transfer coefficient Is seen to rise slightly for about the first 6 min of 
Icing. This Is probably due to the. Increased roughness caused by Ice accre- 
tion. The heat transfer then falls off as the flow Is reduced due to the 
higher pressure drop. Figure 16 shows that the overall heat transfer coeffi- 
cient Is much less sensitive than the pressure drop to Ice accumulation. 


CONCLUSIONS 

Air side pressure drop and overall heat transfer characteristics of a 
heat exchanger proposed for use In the NASA Lewis Altitude Wind Tunnel have 
been tested In the Icing Research Tunnel under dry and Icing conditions. The 
objective of these tests was to determine If the proposed heat exchanger could 
meet the requirements of AWT. The conclusions are: 

1. The proposed heat exchanger exceeds the design heat transfer require- 
ments by 26 percent. 

2. Dry pressure drop Is reasonable at 1420 Pa (5.7 In of water) (heat 
exchanger only - not Including wide angle diffuser ("blister") losses). 

3. Frost forms uniformly on all heat exchanger surfaces, consequently AWT 
cannot be started with moist air. 

4. Ice forms only on fin and tube leading edges that are In direct llne- 
of-slght to oncoming stream. 

5. Ice builds up faster at the colder air temperatures - dimensionless 
pressure drop doubles after 10 min Icing with 15 vm drops at 244 K (-20 °F); 
completely blocked In 20 to 30 min at this condition. 

6. Ice builds up faster with smaller drop sizes probably due to freezing 
of water drops Into Ice crystals before contact with heat exchanger surface. 

7. Surface temperature has little or no effect on Ice accretion rate If 
air temperature sufficiently below freezing. 
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TABLE I. - HEAT EXCHANGER SEGMENT PARAMETERS 


Minimum free flow area, A m 

Total air side heat transfer area, A 0 
Free flow area/frontal area, A m /Af . . 

Total fin area, Afj n 

Fin area/total area, Af^ n /A Q 

Refrigerant side heat transfer area, Ar 

Surface area/volume 

Hydraulic diameter, 4r^ 


. . . 0.419 m 2 (4.506 ft 2 ) 
. . 97.2 m 2 (1046.15 ft 2 ) 

0.509 

. . 91.22 m 2 (981.9 ft 2 ) 

0.939 

... 6.08 m 2 (65.46 ft?) 

340.9 m 2 /m 3 (103.9 ft 2 /ft 3 ) 
. . . 0.549 cm (0.018 ft) 


TABLE II. - REQUIRED HEAT EXCHANGER PERFORMANCE AT VARIOUS CONDITIONS 
(-55 °F FLOODED COIL) 


[Total Heat Load = Q fan + Q ice + Q s team/ heat -] 



Mach 

number 

Altitude, 

ft 

Temperature, 

°F 

Heat load 
xlO , 
Btu/hr 

Mass 

flow, 

lbm/sec 

LMTD, 

°F 

Re 

UAxlO" 6 , 

Btu/hr-°F 

1 

0.20 

762 

439.7 

128.77 

5 484 

46.9 

19 608 

2.7406 


.20 

23 163 

439.7 

23.47 

2 353.7 

40.1 

8 645 

0.5839 


.30 

1 710 

521.8 

32.26 

7 607.6 

119. 

24 780 



.30 

24 366 

439 7 

39.6 

3 368.7 

41.1 

12 328 

.9628 


.40 

3 012 

524.5 

68.29 

9 713.5 

123. 

31 368 

.5528 

HI 

.40 

26 012 

439.7 


4 208.6 

42.4 

15 320 

1.4274 

1 

.52 

5 014 

527.9 

120.08 

11 806 


37 740 


8 

.52 


439.7 

88.19 

4 948.1 

44.2 

17 890 

1.9945 

9 

.60 

14 875 

499.2 

114.89 

9 639 


32 045 

1.1381 

10 

.60 


439.7 

107.7 

5 268.9 

45.5 

18 956 

2.3665 

11 

.70 


471.7 

163.66 

7 665.9 

78.4 

26 140 

2.0873 

12 

.70 

33 215 

439.7 


4 584 

46.9 

19 615 

2.7406 


.80 


456.6 


6 595.5 

65.2 

22 862 

2.5744 

14 

.80 

55 000 

439.9 

102.28 

2 239.3 

57.3 

7 680 

1.7841 

15 

.84 


453.3 

169.26 

6 340.5 

62.4 

22 040 

2.7085 

16 

.84 

55 000 

444.6 

108.64 

2 341.9 

62.6 

7 963 

1.7338 

17 

.92 


455.4 

119.98 

2 564.2 

74.2 

8 569 

1.6160 


m = 0.3048 x ft. 

Watt = 0.2929 x Btu/hr. 
kg/sec = 0.4536 x lbm/sec. 
deq K = (°F + 459.67) 1.8. 
W/ C = 0.5275 x Btu/hr-°F. 




















TABLE III . - HEAT TRANSFER AND PRESSURE DROP DATA 


Rdg . 

number 

Re 

Error , 

X 

u, 0 

Btu/hr-ft 2 -°F 

Error , 

% 

25. 

46830.8 

2.4 

8.78 

4.0 

26. 

46396.7 

2.4 

7.34 

7.6 

27. 

46446.7 

2.4 

8.78 

3.4 

28. 

45065.9 

2.5 

8.89 

3.1 

29. 

44902.8 

2.5 

9.62 

3.4 

30. 

26764.1 

7.6 

7.20 

7.9 

31. 

28858.8 

6.1 

7.90 

6.6 

32. 

28473.2 

6.3 

8.02 

6.8 

33. 

27746.7 

6.7 

8.02 

7.5 

34. 

27492.3 

6.7 

7.95 

7.4 

35. 

26190.6 

7.4 

7.87 

8.2 

36. 

25911.9 

7.6 

7.81 

8.7 

38. 

23217.6 

9.4 

7.39 

10.9 

39. 

18387.9 

16.6 

5.95 

17.5 

40. 

18467.2 

17.0 

6.05 

18.0 

41. 

17936.9 

17.7 

5.93 

19.3 

42. 

13411.6 

47.9 

4.80 

48.2 

43. 

12161.6 

33.9 

4.40 

34.5 

44. 

11373.5 

39.5 

4.22 

40.0 

45. 

14512.9 

26.4 

5.15 

27.8 

47. 

41331.0 

3.4 

6.65 

4. 1 

48. 

37648.5 

3.9 

6.55 

4.4 

49. 

33508.6 

4.5 

6.44 

6.0 

50. 

27211.9 

6.7 

5.76 

10.5 

51. 

18555.9 

15.8 

4.72 

16.5 

52. 

10843.6 

45.0 

3.49 

45.3 

54. 

62055.3 

2.0 

11.23 

7.5 

55. 

59257.7 

2.1 

11.09 

7.7 

56. 

57421.9 

2.2 

11.42 

8.3 

57. 

44716.8 

3.7 

11.56 

9.4 

58. 

14054.7 

42.5 

5.80 

44.2 

60. 

55720.9 

2.2 

9.40 

17.0 

61. 

49917.5 

2.4 

7.28 

12.8 

62. 

49804.2 

2.5 

7.83 

10.2 

63. 

45753.4 

2.9 

8.16 

10.5 

64. 

50482.0 

2.4 

8.90 

12.4 

65. 

49132.3 

2.5 

9.24 

11.6 

66. 

45104.4 

2.9 

9.15 

11.9 

67. 

33940.3 

5.1 

6.77 

14.1 

70. 

57796.3 

2.1 

10.91 

8.0 

71. 

58356.2 

2.6 

11.29 

8.7 

72. 

58516.9 

2.1 

11.81 

8.6 

73. 

53722.9 

2.5 

11.47 

8.7 

74. 

45618.6 

3.4 

10.64 

8.8 

75. 

24321.9 

12.8 

7.80 

15.5 

77. 

51376.8 

2.1 

10.95 

8.1 

78. 

56034.0 

2.1 

9.56 

7.8 

79. 

53951.4 

2.3 

9.21 

9.5 

80. 

52041.0 

2.4 j 

9.54 

12.6 

81. 

49378.7 

4.8 

9.40 

10.3 

82. 

51334.1 

2.5 

10.01 

9.4 

83. 

30504.2 

7.0 

7.82 

14.4 

93. 

58491.0 

2.2 

11.31 

8.7 

94. 

57782.4 

2.2 

10.72 

9.0 

95. 

58067.2 

2.2 

11.05 

9.1 

96. 

57032.0 

2.4 

11.39 

10.2 

97. 

55285.0 

2.4 

11.23 

9.8 

98. 

54270.7 

2.6 

11.53 

10.0 

99. 

49907.1 

2.8 

11.51 

10.0 

100. 

39607.4 

4.7 

11.39 

11.5 

111. 

39369.7 

5.0 

9.36 ! 

8.2 

112. 

29691.4 

10.0 

8.58 

12.2 

113. 

42184.3 

4.3 

9.46 j 

8.0 

114. 

50626.1 


10.27 

7.6 

115. 

55273.3 

2.6 

10.59 

7.5 

116. 

58635.9 

2.4 

11.03 

7.6 

117. 

32657.1 

7.4 

9.46 

10.3 

118. 

49622.6 

3.2 

10.59 

7.8 

119. 

60176.5 

2.2 

11.12 

7.6 

W / m 2 °C 

= 0.1761 

x Btu / hr - ft 2 -° F . 



AP S 

q 

Error , 


Error , 

Mce* 

min 

tfrost 

min 


0.4 

42.17 

1.4 


11. 

0 

■ It 


42.37 

2.0 



16.0 

■ 


40.22 

1.6 



21. 

0 

14.48 


39.25 

1.7 



26.0 

14.73 


38.38 

1.8 



31. 

0 

14.89 

1.3 

38.83 

1.5 



6.0 

14.16 

0.9 

38.60 

1.9 




0 

15.77 

0.9 

37.86 

2.1 



16.0 

17.91 

0.8 

37.19 

2.5 



22.0 

19.30 

0.8 

36.93 

2.5 



27. 

0 

22.46 

0.7 

36.00 

2.6 



32.0 

24.07 

0.7 

35.59 

2.9 



41. 

0 

12.79 

1.6 

37.20 

2.8 



7. 

0 

24.27 

1.3 

36.98 

3.0 


1 

12.0 

27.02 

1.2 

36.48 

3.2 



17. 

0 

31.47 

1.1 

36.19 

4.0 


. 1 i 

22.0 

61.53 

1.0 

35.27 

3.4 



27. 

0 

80.15 

0.9 

34.91 

3.8 



32.0 

97.64 

0.9 

34.69 

4.0 


■ 

37.0 

63.26 

0.8 

35.42 

4.9 


■ 

42.0 

18.95 

0.4 

41.79 

1.4 


1 

5.0 

26.16 

0.4 

43.21 

1.4 


1 

10.0 

37.12 

0.3 

43.70 

1.6 



15.0 

58.81 

0.6 

43.87 

2.6 



20.0 

125.63 

0.3 

43.09 

1.9 



25.0 

377.61 

0.3 

41.91 

1.9 



30.0 

10.66 

0.4 

7.36 

6.1 



0.0 

12.92 

0.3 

7.28 

6.4 




15.17 

0.3 

6.45 

6.9 

10.0 



34.07 

0.3 

5.86 

7.9 

20.0 



345.78 

0.2 

5.50 

9.0 

30.0 



10.18 

0.5 

5.13 

7.6 

0.0 



14.55 

0.4 

4.35 

7.5 

5 

0 



15.48 

0.4 

4.49 

7.4 

10.0 



19.21 

0.3 

4.23 

7.9 

15.0 



15.06 

0.3 

3.46 

9.5 

20.0 



16.44 

0.3 

3.65 

9.1 

25.0 



21.64 

0.3 

3.56 

9.6 

35.0 



42.98 

0.3 

3.22 

10.8 

50.0 



11.73 

0.4 

7.03 

6.6 

0.0 



11.86 

0.4 

6.51 

7.1 

2.0 



12.81 

0.4 

6.65 

7.3 

5 

0 



16.95 

0.3 

6.60 

7.4 

10.0 



26.88 

0.3 

6.67 

7.3 

15.0 



105.20 

0.3 

5.97 


20.0 



10.74 

0.5 

7.33 

5.0 

0.0 

6.0 

11.21 

0.4 

6.82 

6.1 

0.0 

0.0 

13.04 

0.4 

5.38 

7.1 

5.0 



14.95 

0.3 

4.77 

8.3 

10.0 



17.46 

0.4 

5.28 

7.7 

15.0 



16.83 

0.3 

5.19 

7.6 

20.0 



60.04 

0.6 

4.38 

10.6 

30.0 



11.80 

0.4 

6.40 

7.1 




12.71 

0.4 

6.28 

7.2 




13.14 

0.3 

5.90 

7.5 

10.0 



14.52 

0.4 

5.78 

8.1 

,s 

.0 



15.96 

0.3 

5.76 

7.9 

20.0 



17.07 

0.5 

5.32 

8.4 

25.0 



22.60 

0.3 

5.43 

8.6 

35.0 



42.74 

0.3 

5.08 

9.8 

45.0 



13.36 

0.8 

9.09 

5.8 




15.11 

1.3 

8.87 

6.1 


m. 



16.04 

0.6 

8.78 

5.9 

■ 

I 



15.28 

0.4 

8.84 

5.9 

■ 




15.02 

0.4 

8.66 

5.9 


■ 



15.24 

0.4 

8.55 

6.2 





13.31 

1.3 

8.38 

6.3 





14.75 

0.5 

8.49 

6.2 



1 


15.09 

0.4 

8.41 

6.1 

■ 

m 
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Figure Z - AWT icing goal and FAR 25 icing envelope. 



r- Location of heat exchanger 
/ Module test duct 
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Figure 3. - Schematic of icing research tunnel. 



Figure 4. - IRT refrigeration system with AWT heat exchanger segment. 





(a) AWT heat exchanger test duct. 
Figure 5. 
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(b) Duct inlet, heat exchanger and upstream rake. 
Figure 5 . 
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Figure 8. - Totai and static pressure drop-dynamic head 
ratio. 



Figure % » Heat exchanger overall heat transfer coefficient 
based on air side area. 
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Ice accretion time, min 

Figure 15. - Effect of heat exchang- 
er surface temperature on pres- 
sure drop due to ice accretion 
for 244 K (-20 °F), 15 pm drops, 
1.49gm/m 3 LWC (AWT equiv- 
a lent). 



Ice accretion time, min 


Figure 16. - Effect of ice accretion 
on heat exchanger overall heat 
transfer coefficient for 294 K 
(-20 °E>, 15 pm drops and 1.49 
gm/rrr LWC (AWT equivalent). 
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